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Abstract 

Understanding the rich dynamics of open quantum systems is of fundamental interest to quantum 
O • control and quantum information processing. By considering an open system of many identical 

two-level atoms interacting with a common bath, we show that effects of system-bath correlations 
are amplified in a many-body system via the generation of a short time scale inversely proportional 
to the number of atoms. Effects of system-bath correlations are therefore considerable even when 
each individual atom interacts with the bath weakly. We further show that correlation-induced 
dynamical effects may still be suppressed via the dynamical decoupling approach, but they present 
a challenge for quantum state protection as the number of atom increases. 
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I. INTRODUCTION 

The rich dynamics of open quantum systems continues to attract great interests [lj • Two 
main motivations for now are to achieve better quantum control and to better understand 
possible quantum effects in large systems. However, the general complexity of open quantum 
systems puts exact solutions out of reach. A variety of approximations or assumptions are 
therefore needed, with their validity under close scrutiny in recent years due to fascinating 
experimental advances in, e.g., cold-atom physics and photonics. One common assumption 
in treating open quantum systems, often referred to as "factorized initial conditions", is 
that the system and the environment are initially uncorrelated. This can be justified for 
weak system-bath coupling because the system has a negligible impact on the bath statistics 
On the other hand, for moderate and strong system-bath coupling, which is the case 
in some realistic situations of great experimental interest (e.g., light-harvesting systems, 
super-conducting qubits, and atom-cavity systems), initial system-bath correlations (SBCs) 
should be accounted for 0-7]. 

Pioneering studies of SBCs in single-body systems, e.g., a sing! 
oscillator in a thermal bath, have been fruitful js-23| (Refs. js, 24] are notable exceptions 
involving two spins). The main purpose of this work is to extend the investigation to a 
system of many identical two-level atoms, where novel collective phenomena might occur 



e spin or a single harmonic 



(one example is super-radiance 25|). Specifically, we consider a system of many identical 
two-level atoms, each atom weakly interacting with a common bath. Using an exactly 
solvable model, we show below that effects of SBCs may dramatically increase with the 
number of atoms N, leading to a previously unknown time scale t c inversely proportional 
to N. For a large system with many atoms, t c becomes very small and SBCs manifest 
themselves by inducing rapid oscillations in physical observables. Three implications of our 
findings are in order. First, the collective dynamics of many two-level atoms can be strongly 
affected by SBCs even when each individual two-level atom interacts weakly with the bath. 
Second, this fact may be exploited to amplify and gauge SBCs. Third, a small t c presents 
challenges for quantum state protection via dynamical decoupling (DD) techniques. In 
general, it is found that SBCs force us to apply more frequent and more efficient DD pulses 
to freeze the quantum evolution. 

In Sec. II, we present our theoretical findings of SBC effects in a model system consisting 



of N two-level atoms in a common bath. Section III presents computational examples. We 
then discuss in Sec. IV the implications for the dynamical decoupling technique and how it 
may be used to suppress SBC effects. Section V gives a brief summary. Readers interested 
in the technical details of our derivations should refer to Appendices A-D. 



II. SYSTEM-BATH CORRELATION EFFECTS IN A SYSTEM OF N TWO- 
LEVEL ATOMS IN A COMMON BATH 

A collection of N identical two-level atoms interacting with a common bosonic bath may 
be described by an extended spin-boson Hamiltonian H total = H s + H B + H mt (setting h = 1 
throughout) 3,0, with 

H s = u J z + 6J X + X Jh (!) 
H B = J2^kblh, (2) 

H int = 2J z J2(9th + 9kbl). (3) 

k 

Here J x , y , z are the collective spin operators with J x + Jy + J z = y(y + 1), is the 
energy bias, x the interaction between the atoms, 5 is the tunneling amplitude, and H B is a 
collection of boson modes or harmonic oscillators (with zero-point energy dropped). Such a 
total Hamiltonian H tot3i \ can model a two-mode BEC interacting via collisions with thermal 



atoms or phonon excitations [28j. We stress that, other than the self-interaction term x-J^ 



29| to a 



H tota i is nothing but a straightforward extension of the standard spin-boson model 
model of many spins interacting with a common bath. 

Rather than switching on the system-bath interaction H int at a particular instant, in most 
physical situations the system and the bath have interacted for a long time beforehand. Our 
starting point is then a thermal equilibrium state for the system and the bath as a whole 
at temperature T, i.e., p oc exp(— [3H total ) with = Noticing the energy contribution 
by H- mt , one finds that the associated reduced state of the system (bath) is not really given 
by p^ q oc e~^ Hs (p c B oc e _/3i ^ s ) {?], S]- Instead, states of the system and of the bath are 
correlated for a non- vanishing H [nt . Now, if at time t = 0, the system is prepared in a pure 
state |^) via a projective measurement, then the initial state of the system and the bath as 



a whole is given by 



, (ih\ e~P H «*<* lib) , 
p(0) = |V) M ® — ^ — , (4) 



where Z is a normalization factor. This initial state should be compared with the usual 
uncorrelated (unphysical) initial state, 

Evidently, for a nonzero i?i n t> we have p(0) 7^ p dir (0). That is, the physical initial state 
in Eq. (j4j) has correctly accounted for the system-bath interaction during the past. Conse- 
quently the initial bath state p_e(0) = \ e~^ Htotal /Z depends on H int as well as the state 
preparation of the system and is therefore not a canonical equilibrium state for the bath. 
Previously, how the difference between p(0) and p dir (0) impacts on the ensuing dynamics was 



studied for a damped harmonic oscillator (see, e.g., Ref. [13]) and a single two-level system 



undergoing pure dephasing [22[. In Appendices C and D, we also consider initial states of 
the form Qe~ l3Htota - 1 QyZ, where Q is a unitary operator acting on the system Hilbert space, 
and we argue that this initial state for large N leads to similar results as p(0) does. 

In general, the dynamics of H tota \ starting from p(0) or p dir (0) does not have analytical 
solutions. To obtain analytical solutions from which important insights may be gained, we 
set the tunneling parameter 5 in H$ to zero, yielding a pure- dephasing problem for collective 
spin states. In the context of a two-mode BEC, such a situation arises if intermode coherent 
tunneling is made to vanish and if intermode mixing collisions are negligible. Purely for 
convenience, we shall assume \ = 0, which may be achieved via Feshbach resonance. Up to 
a unitary transformation, one can obtain equivalent situations if = and if inter-mode 



mixing collisions dominate 
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Working in the basis of J z eigenstates, denoted by \m) with J z \m) = m \m), we first find 
the system's reduced density matrix for the uncorrelated initial state p dir (0) (see Appendices 
for detailed calculations), 

\ps(t)]mn = [ps(0)Ue-^ m -^e-* A «( m2 -^ 

x e -^)^?\ (6) 

where the decoherence factor is given by 



7M = £4| 9t P<I^«coth(^), (7) 



and 



A(i) 



sin(u k t) - u k t 



(8) 



The factor exp[— iA(t)(m 2 — n 2 )t] arises because the common bosonic bath assists in gen- 
erating an indirect atom-atom interaction. Note that at sufficiently short times for which 
sin(o;fct) ~ cu^t, we have A(t) « 0. 

By contrast, for the physical initial state p(0) in Eq. (J4]), we find 



[ps(t)} mn = [ps(0)\ 

x e - 7(t)(m - n)2 *F^ n (t) 



g— iu>o(rn— n)t — iA(t)(m 2 — n 2 )t 



where 



with the probability 



and 



= 2> e ' 



-i[2/(n-m) £ fe 4| 9fc | 2 sin(u; fc *)/u;g] 



2 e -f3LO l e l3l 2 C 



Ez(l(W| 2 e 



2 p -Puiol e l3l 2 Cj ' 



C = y^4|^ fc | 2 /a; fc . 



(9) 



(10) 



(11) 



(12) 



Comparing Eq. (jSJ) with Eq. (jUJ), it is seen that effects of the SBC on the dynamics are entirely 
captured by the factor F^ n (t). To understand this, we first examine the non-canonical bath 
state ps(0), i.e., the reduced bath state at time zero. Up to a normalization factor, we 
obtain 

p B (0)o<J2^ l \m\ 2 e-^\ (13) 



where 



H® = H B + 2lJ2(9*kh + 9kbl: 



(14) 



_Ott0) 

Two observations can be made here. First, the bath is prepared in the state e~> B with a 
probability related to the system projection amplitude (l\ip)- Second, H^' can be interpreted 
as a collection of harmonic oscillators, each of which is under a 'force' proportional to 21. Due 
to this force exerted by the system, the actual equilibrium position of the bath oscillators 
will be displaced. To make this clearer, we define displaced harmonic oscillator modes 



Bk,i = bk + 2lg k /uk, 



(15) 
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from which we have 



(16) 



The initial bath state is seen to be a mixture of different components: each component 
is a canonical equilibrium state oc e -/3 Efc aJfe ' B fc,i B M for a collection of harmonic oscillators 
displaced by 2lgk/uik, with the probability pi defined in Eq. ffTT]) . 

An interesting physical picture then arises. At time zero a projection of the system 
on state breaks the equilibrium state of the system and the bath as a whole and the 
bath is instead prepared in a mixture of many components ~ e -PT.kUkB kl B k ^ ^ w ^ 
component assuming a certain J z eigenstate for the system. For the |m) component of 
the component e _/3 ^'= Wfc ' B fe. I ' Bfe i of Pb(0) finds its equilibrium condition no longer satisfied for 
I ^ m. As a result the collection of displaced harmonic modes start to oscillate around their 
new equilibrium positions defined by 2mgi s ./uJk- For another \n) component of \ip), the same 
mechanism works but with a different degree due to new equilibrium positions at 2ngk/uJk- 
It is such type of bath motion that yields the correction factor F^ n (t) in Eq. fl9]). 

Because e^ [2c changes rapidly with / G [— N/2, N/2], for large TV only the terms with 
probabilities p ± N may contribute to F^ n (t). Further, for NftujQ ^> 1, we have e NI3ul0 ^ 2 ^> 
e -N/3ui /2 an( j h ence p_ jv ^> pN_ for a generic state \ip). This yields p_ n « 1 in our many-body 
system, thus reducing F^ n (t) to 



In particular, at sufficiently short times, sin(w fc t) « u k t and then F^ n (t) « e l [ N{ - n ~ m ^ t ^k i \9k\ 2 /^ k ] 
That is, F^ n represents a phase factor building up with time at a rate of N(n—m) 4\gh\ 2 /uik- 
Note that during the same time window, A(t) originating from bath-assisted atom- atom 
interaction still stays close to zero due to sm(ukt) ~ u^t. For physical observables not 
diagonal in the \m) representation, the time dependence of F^ n leads to a characteristic 
time scale (after setting \m — n\ = 1): 



Remarkably, t c is inversely proportional to the number of two-level atoms. As an example 
we consider an Ohmic spectral density for the bath, with 



mn\ J 



C (t) RS e i[N(n-m)J2 k Mgk\ 2 sm(u> k t)/u)l] 



(17) 




(18) 




(19) 
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and 

J{u) = Goue'^. (20) 

Then one finds t c = (NGuc)' 1 , a time scale determined by parameters from both the system 
(N) and the bath (G and u c ). On the other hand, at later times, the time dependence of 
F^ n (t) weakens as the oscillations of the bath oscillators around their new equilibrium 
positions start to dephase. For the Ohmic spectrum, F^ n (t) at long times is found to 
approach e l ^ NG ^ n ~ m ^ . Analogous results are found for other spectral density functions 
J(oj) with an exponential cutoff. 



III. NUMERICAL EXAMPLES OF AMPLIFIED SYSTEM-BATH CORRELA- 
TION EFFECTS 

In this section we turn to a concrete example for which we investigate the dynamics of a 
scaled observable j x = 2(J X )/N. At time zero, the system is projected onto the initial spin 
state \ip) = e~ l ^ Jy \N/2). Using the uncorrelated initial state in Eq. (jSJ), we obtain 

j x = e" 7(f)i Re [e l[wo+A(t)lt X(t)] , (21) 

1 
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FIG. 1. (color online) j x vs t without a bath (upper dotted, blue), with an Ohmic bath but using a 
factorized initial state (dot-dashed, magenta), and with an Ohmic bath and including SBC effects 
(solid, black). N = 2000, ujq = 0.1, G = 0.001, = 1000, and uj c = 10 Q- Inset shows the parallel 
results if N = 1, where three lines become almost indistinguishable. 
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FIG. 2. (color online) Same as in Fig. [H but now with N = 20000. The rapid oscillations in 
j x clearly demonstrate the amplification of SBC effects achieved by an increase in the number of 
particles. 



with 



w at or 



2imtA(t) 

N^2 N 

m 



.f + m )(f + Ll)(f- m )(^ + m + 1 )' <22) 

Details can be found in Appendix D. By contrast, the true dynamics with the physical initial 
state in Eq. (]1J) gives 

Jx = e-rWtRe [e^+^X^F^t)] , (23) 

with F£(t) ~ ^NGuct snor t times for an Ohmic spectrum defined above. In Fig. [TJ we plot 
the time dependence of j x for N = 2000, in the absence or presence of a bath with G ^ 1 



(so that each individual two- level atom interacts with the bath very weakly |26[). For the 
shown time period in Fig. [H j x in the absence of the bath hardly changes due to a finite 
ojq. In the presence of the bath but without including SBCs, j x stays positive but decreases 
at a faster rate due to the bath-assisted atom-atom interaction. With SBCs accounted for, 
completely different qualitative behavior of j x is observed: it rapidly changes from 1 to 
— 1 and then gradually returns to a value close to zero. Note that, as analyzed in theory, 
such rapid change in j x occurs before atom-atom indirect interaction (as captured by A(t)) 
takes effect. The inset also shows the parallel results if N = 1, with all other parameters 
unchanged. As expected from G 1, in that case all lines are almost on top of each other 
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and hence no SBC effect can be seen. Because the decoherence function j(t) is independent 
of N, the results in the inset also hint that decoherence for the N = 2000 case is insignificant 
for the shown time scale. Therefore, SBCs impact on the dynamics long before the onset of 
decoherence. 

To emphasize the role of N, in Fig. [2] we plot parallel results for iV = 20000. There SBCs 
induce more drastic oscillations in j x , followed by slower oscillations as the time dependence 
of F^ n (t) weakens. Clearly then, at short times an increasing number of atoms enhances 
the oscillation frequency in physical observables, thus amplifying the SBC effects. 



IV. SUPPRESSION OF SYSTEM-BATH CORRELATION EFFECTS BY DY- 
NAMICAL DECOUPLING 

In this section we discuss the implications of this work for DD, which is one main ap- 



proach to quantum state protection [31H33I] . DD in single-spin systems has found enormous 
applications. It is therefore highly desirable to extend DD to systems describable by col- 
lective (large) spins J x , y ,z |3j 



351 ] . We return to our pure-dephasing model introduced in 
Sec. II. Because the details of a bath spectrum is unknown in general, the precise form of the 
correction factor F^ n (t) induced by SBCs or of the A(t)-related phase factor due to bath- 
assisted atom-atom interaction is unavailable in general. It is tempting to wait for F!^ n {t) 
to saturate such that its time dependence is out of the picture. However, as indicated by 
our theory and shown in Fig. [2J before reaching that regime the bath assisted atom-atom 
interaction would have already changed the state. So in order to protect or store a given 
many-body state, the time dependence of F^ n (t) and A(t) must be suppressed. In general 
the task of DD becomes three-fold: to suppress the decoherence factor 7(t), the A(t)-related 
phase factor, and the correction factor F^ n (t). For large N, it is found that DD control 
pulses need to first compete with the correlated-induced time scale t c oc 1/N. 

Consider now Nd instantaneous 7r-pulses of J x applied to our system at times ti, with 
1 < I < Nd- Upon application of one such pulse, one has, in the frame of the applied pulses, 
J z — > — J z . It is hence convenient to introduce the so-called switching function fit), with 

N d +l 

f(t)= ^(-l^t-fj.Oflft-t), (24) 
i=i 

where 9{t) is the Heaviside function. With the assistance of f(t), one can still find the time 



dependence of [ps{t)]mn or of an observable such as j x analytically, using the same technique 
as used in previous cases without DD pulses. Take [ps(t)]mn under DD as an example. The 
general form of Eq. ([§]) still holds, but now with u (t) changed to u (t), A(t) changed to 
A(t), 7(t) changed to j(t), and F^ n (t) changed to F^ n (t). Specifically, we find 

£ = ^ f f(t')dt', 
1 Jo 

A(t) = V %^ / dt x f 1 dtzfitjfih) sm[cu k (t 2 - t,)}, 
k t Jo Jo 

m ^ Mg ^\!Mi coth ^y (25) 



and 



Here 



F c (t) « e " i [ Ar(n " m) ^ fe4|9fc|2lm[/(wfc ' <)1/tJ fe]. (26) 

N d 

f(co k , t) = l + (-l) N * +l e luJkt + 2 ^(-l)' e ^ (27) 

i=i 

is determined by f(t) via the following relation 

f(u k ,t) = -iu k f dt'e^'fit'). (28) 
Jo 

With these explicit results, effects of DD pulses on the dynamics with SBC effects can be 
examined in detail. Of particular interest is F^ n (t) because it is indicative of the impact of 
an applied DD sequence on SBC effects at early times. 

In Fig. [3] we illustrate the effect of applying DD pulses for N = 20000 for a final time 
t — 0.1. We first investigate the usefulness of equidistance control pulses (bang-bang control) 
with a pulse interval r 3l| • Note first that if we neglect SBCs, then there would be no need to 
apply any DD pulses because the evolution of j x is negligible for the considered period. With 
SBCs accounted for, an application of Nd = 4 pulses with r = 0.02 is found to be insufficient 
for state protection. Indeed, for the system parameters used we find t c = (GNuc)" 1 = 0.005, 
which is far smaller than r = 0.02. To compete with t c , we then use r = 0.002 < t c , for 
which the evolution is successfully frozen (see Fig. [3]). 

To avoid the need of a high pulse repetition rate, we next consider the celebrated, more 



effective, Uhrig's dynamical decoupling (UDD) sequence [32|,|33| with unequal pulse intervals, 
where the pulse timings are given by tj = tsin 2 ^ 2J ^ 7r +2 ) ■ Dramatically, via only Nd — 4 
UDD pulses, the j x value at the final time t — 0.1 is already recovered to its initial value. 
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FIG. 3. (color online) j x vs time with SBC effects and under bang-bang control pulses, for r = 0.02 
(solid, black) or r = 0.002 (dotted, black). The parallel results without any control pulses (dotted, 
magenta) and with four UDD pulses (solid, magenta) are also shown. Parameters are the same as 
in Fig. [21 with the final time chosen as t = 0.1. 



A careful look into the above expression for F!^ n (t) explains why this is so. The original 
motivation of a UDD sequence is to suppress the decoherence function 7(t) by minimizing 
\f(cuk,t)\ to its A^-th order in time. So by construction, the time evolution of F^ n (t) is 
also optimally suppressed by UDD because a minimized \f(uk,t)\ automatically yields a 
minimized lm[f(uk,t)} that enters into F^ n (t) in Eq. (1261) . That is, here the correction 
factor F^it) is already optimally suppressed for an unknown spectrum with an exponential 
cutoff. Additional computational studies indicate that the A(t) term may be also well 
suppressed by DD pulses, though this is not of interest here as our main concern is to 
suppress SBC before bath-assisted atom-atom interaction has any considerable effect on the 
dynamics. Our conclusions are as follows. In locking a many-body state in our model here, 
the previously unknown time scale t c induced by SBC calls for the use of DD control pulses 
long before decoherence and bath-assisted atom-atom interaction becomes important, with 
a UDD sequence found to be an optimized choice (assuming the detailed form of the bath 
spectrum is not available). 
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V. SUMMARY 



In summary, because SBC effects are amplified by the number of particles in a many-body 
system, they can be important even when each individual particle interacts with the bath 
weakly. Interestingly, though SBC effects originate from system-bath interaction in the long 
past, they may still be suppressed by dynamical decoupling so that a prepared many-body 
state is well protected. Nevertheless, reaching this goal calls for more effective control pulses 
applied within a shorter time scale. Our results should also be of interest to other subtopics 
in open many-body systems by considering SBC effects neglected before. 
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Appendix A: Exact unitary evolution operator 



The dynamics for our model described by Eqs. ([I])-® can be exactly solved if 5 = 0. 
Purely for convenience we also assume x = 0- We first transform to the interaction picture. 
In this picture, the Hamiltonian becomes 



Hj(t) = e Ws+H B )t Hinte -i(H s+ H B )t^ 



(Al) 



Similar to the treatment in Ref. 



we next find the time evolution operator Ui(t) corre- 



sponding to Hj(t). The Magnus expansion 



tells us that 



Uiit) = exp 



1=1 



(A2) 



where 



A 1 (t) = -i [ dt x Hi{h), 
Jo 

* Jo Jo 



(A3) 
(A4) 
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It is straightforward to find 

Ai = -z / dtiHjfa), 
Jo 

= J z J2(bla k (t)-b k al(t)), (A5) 



where 



Further using 



we obtain 



where 



, , 2g k (l - e tuJkt ) , K . 

a k (t) = '-. (A6) 

[if/fa), flj(t 2 )] = -8tJ 2 \9k\ 2 MMh - h)}, (A7) 

k 

A 2 = -[ dti I dt 2 SiJ 2 y2\g k \ 2 sm[u k (t 1 -t 2 )], 
z Jo Jo k 

= -iJ 2 tA(t), (A8) 

A(t)E^%l' [afaM r^ - (A9) 

Since this is a c-number, the higher order terms in the Magnus expansion are all zero. The 
exact unitary time evolution operator is hence found, i.e., 

U(t) = e- Mt e- iHBt U!{t), (A10) 

where 

Z7/(t) = exp{ J z Y} h W{t) - b k a* k (t)} - U 2 tA(t)}. (All) 

k 

For later calculations, let us first consider the reduced density operator of the system 

p g (t) = Tr B [U(t)p(0)tf(t)}, (A12) 

where p(0) is the density operator of the system and the bath as a whole. We find it useful 
to write the reduced density operator of the system in terms of the standard J z basis as 

[ps(t)] mn = Tr s> B[U(t)p{0)U\t)P nm }. (A13) 

Here P nm = \n) (m\, \n) being the eigenstate of J z with eigenvalue n. Introducing the 
Heisenberg picture operator P nm (t) as W (t)P nm U(t), we have 

[ps(t)]mn = Tr s>B [P nm (t)p(0)}. (A14) 
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Calculations of the explicit form of P nm (t) are straightforward because the time evolution 
operator U (t) is already found. In particular, we obtain 

P nm (t) — ^^oKrn-n) e -iA(t)t(rn 2 -n 2 ) e -R nm {t) (A15) 

where 

R nm {t) = {n-m) Yy>k a k{t) - b k a* k {t)}. (A16) 
k 

It follows that 

[ps(t)]mn = e-^*( m -^e- jA W'( m2 - n2 )Tr 5iB [e" if -«P n , m p(0)]. (A17) 

This is a general result because it applies to an arbitrary initial density p(0) (the unphysical 
state p dir (0) or the physical state with SBC accounted for). 

Appendix B: Dynamics with uncorrelated initial states 

Here we consider unphysical decorrelated initial states, i.e., 

p dir (0)=p 5 (0)®p B , (Bl) 
where p# = e with Zg = Trg[e~P HB ]. Then, 

[ps(t)}mn = [p5(0)] m „e-^( m -^e- iA W( m2 -" 2 ) t Tr B [e- fl ™WpB]. (B2) 

We now simplify Tr #[e~ - Rnm W p B ] = ^-^C)^ where the average is taken with respect to the 
thermal bath state at equilibirum. Although this is a standard result, for self-completeness, 
we explain the steps in detail. Since the modes are independent of each other, we can write 

^g--Rnm(t)^ _ J~J^ e -(n-m)[&£a fc (t)-6 fc a£(t)]^ (33) 
k 

For an operator A which is a linear combination of creation and annihilation operators, we 
have (e A ) = e^ 2 ^ 2 . Using this identity, we have 



/g — -Rnm (t) \ 



) = JJexp -\{n - mf\a k {t)\ 2 {b[h + b k b\) 



Jexp 



2 

~{n-m) 2 \a k {t)\ 2 (2n k + l) 



(B4) 
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Further using the definition of ak(t) and the Bose-Einstein distribution, we find that 



Tre[e Rnm{t) p B ] = exp 
This then yields 



^(m - n)' 2 A\g k 

k 



2 [1 -cos(u; fe t)] ^ ( /3uj k 



coth 



cur 



[Ps{t)}mn = [Ps(0)] mn e 



-iuio(m—n)t^—iA(t)(m 2 —n 2 )t^—-y(t)(m—n) 2 t 



(B5) 



with 



7 



2 (1 - cos(a; fc t)) /3w fc 

9k\ o cotn 

out 



(B6) 



The factor e~ 1<yt ^ m ~ n ^ 2t describes decoherence and the factor e -* A (*)( m2 -™ 2 )* describes the 
indirect atom-atom interaction induced by the common bath. 



Appendix C: Dynamics with system-bath-correlated initial states 



Here we consider (physical) correlated initial states of the form 

p(o) = ^^n r e- pH ™*nt, 

r 

with 



Tr 



S,B 



(CI) 



(C2) 



22| from a single 



The Q r operators above are assumed to be acting on the system only and their explicit forms 
will be specified later. To solve for the dynamics starting from such an initial state, we may 
use a polaron transformation technique, generalizing the results of Ref. 
spin to many spins, or we may use displaced harmonic oscillator modes. The latter method 
is used below since it is physically more transparent. 

We first simplify Z. By introducing a completeness relation, we find that 



r I 



EE 

r I 



'i\ntn r \i)Tr B [e- pH ^] 



where we have defined 



(C3) 



(C4) 
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Using the displaced harmonic oscillator modes, 

B Kl = h + ^, (C5) 
Bl = b\ + (C6) 

we find 

Z = E (^l^rlO e^Zs, (C7) 

r J 

where C = E fc ^^f > and in this case Z B = Ti B [e^^ hB ^ Bk ' 1 }. 

We now substitute Eq. ([Cl]) in Eq. flATTl) and again introduce a completeness relation. 
To proceeed, Ti B [e~ Rnm ^e~ )3H B > '\ needs to be simplified. As before, this can be done using 
displaced oscillator modes. It is straightforward to show that 

Rnm(t) = (n-m) 5>*(*)4,i - <x*k(f) B k,i] + i&Kt), (C8) 
k 

where 

$(i = 2(n-m)/$(t), (C9) 
$(t) = E%^ sin M)- ( C1 °) 

We then find that 
which leads to 



iojo(m— n)i — iA(t)(m 2 — n 2 )t — j(t)(m— n) 2 t 



[ps(t)}mn = e-^ m - n ^e- l ^ m ~ n ]t e 

J2 r E* ((WiP nm n r \i) e -^Ut) e -^ i e ^c 

x \ L (C12) 

E r E* ((^^r|i}e-^e^ 2c ) 



Noting that 



[Ps(0)}mn = ^7— , ^, (C13) 



one finally arrives at 

[Ps(t)] mn = [psmmne-^- n)t e -iA(t)( m °-n*)t e -y(t)(m-n)H 



E r E* </|^Pn m ^|/> e-^We-^'e^ c 
x ^7 v s-^-- (C14) 
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1. State preparation via projective measurement 



In this case, 

p(0) = Ip^e-^otaip^ ( C15 ) 
with P.jp = Using our general results above we obtain 

[ps(t)} m n = [p^Ol^e-^^-^'e-^W^ 2 -" 2 )^-^^™-™) 2 * 

X Ei(l(^)l 2 e-^°'e^ 2c ) ' (C16) 

This can be rewritten as 

[p S {t)]mn = [p5(0)] m „e-^(— n^'^M^i), ( C 17) 

where F mn {t) is already given in the main text. 
2. State preparation via unitary operations 

Alternatively, instead of performing a projective measurement on the system, we may 
first cool the system and the bath to a desired low temperature. We then perform a unitary 
operation on the system to approximately arrive at some initial state. In general, initial 
state of the system and the bath prepared in this manner can be written as 

p(0) = lfi e - /3Htotal fi t , (C18) 

Zj 

with Q being a unitary preparation operator acting on the system Hilbert space. We then 
have 



[ps(t)]mn = ^(OJl^e-^^^e-^W^-^e^Wf"-")^ 
J2i [(Wn) (m\n\l) e -i^Ut) e -^ i e pi*c 
X (m\n\l) e-^i e ^ 2 c\ • 



This can be rewritten as 

[ps(t)]mn = [p5(0)] m „e-^( m -^ t e- lA W( m2 -" 2 ) i e-^ t )( m -") 2t 

with 



(I) p-i^JLit) 
nm 



(C20) 



{l\Q)\n) (m 









(C21) 

Note that finln need not be real here (nevertheless, by its definition only, we still have 

y~) l pnrn = !)• 
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Appendix D: Calculating observables 



To evaluate physical observables using the system's reduced density operator, let us first 
evaluate the functions j(t), C and A(t). As usual, we take the continuum limit of the 
bath modes, whereby the summations over the bath modes are replaced by integrals via the 
rule 

POO 

J2M9k\ 2 C(u k )= / duJ(u)C(u), (Dl) 
k Jo 

with J{uj) being the spectral density For this work, we choose Ohmic spectral density, that 
is. 

J{uj) = Gue^^. (D2) 



We can then evaluate the integrals as 



22 



and 



C = Gcu c , (D3) 

$(t) = Gtan" 1 ^), (D4) 

7(*)=7vac(*)+7th(*)> ( D5 ) 

7 V ac(t) = ^ln(l + ^ 2 t 2 ), (D6) 

7^)^[lnr(l + ^)-iln|r(l + A- + |)H, (D7) 

m = - t m)-ct). (D8) 



Expressions for sub-Ohmic and super-Ohmic spectral densities can also be found in Ref. 22 1 
that treated a single two-level system in a bath. 

To evaluate the expectation value of J x , we first calculate (J + (t)), where J + = J x + iJ y 
is the standard angular momentum ladder operator. Note first that 



(J+(t)) = ^Pmn(t)(J+W (D9) 



Using 



one obtains 



(n\ J+\m) = \ \ ( — - m J ( — - + m+ 1 ] 5„, m +i- (DIO) 



(J+(t)) =J2P^+i(t)J -mj (y + m + lj. (Dll) 
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We then have 



where 



jx = I (J+(t)) = ^MWXW), (D12) 



For the initial state = e~ l ^ Jy |^) prepared by a projective measurement, it is an 
eigenstate of J x with eigenvalue y. Then, 

[ps(0)] m , m+ i = (m\e~^ J y\N/2) (N/2\e i * J "\m + l) . (D15) 

Note that ^m|e -iaJy |m') is a Wigner "d-matrix" element. Using the Wigner formula, we 
obtain 



Then 



W0)U, = ^( f ^)( f / m + 1 ) ( (D17) 

from which we have 

'W-^?^— ^L)(* + : +1 )(t-»)(t— 0- (D18) 

which is also given in the main text. The observable j x is then given by 

3x = e- i7(i) Re [e i[uo+ *® ]t X(t)FZ(tj\ . (D19) 

As to F£(t) [see Eq. (ID14j) ] that involves the summation over pi, it can be seen that for 
large TV, F£(t) ~ e iAr *(*)_ That is, due to the exponential factors e~^° l and e l3l2c , only the 
I = —N/2 term makes a dominating contribution if N 3> 1. 

We now comment on what happens if, instead of a projective measurement, we use a 
unitary operation. It can be easily shown that 

[ps(0)W cx ( m M l ) (lP\n) e -^ l e^ c . (D20) 
n 
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Once again the term with / = — y dominates, so that 

[ps(0)U « (m|0| - iV/2) (-JV/2|fit| n ) (D21) 

By setting the unitary operator to be e l % Jy , we see that our initial state is approximately 
the same as in the case of state preparation via projective measurement, with F£(t) ~ e lN ®(i) 
for large N. As such, for large N the j x dynamics are very much the same for the above- 
mentioned two realizations of initial state preparation. 
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